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Higgs portals to pulsar collapse 
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Pulsars apparently missing from the galactic center could have been destroyed by asymmetric 
fermionic dark matter {mx = 1 — 100 GeV) coupled to a light scalar (m^ = 5 — 20 MeV), which 
mixes with the Higgs boson. We point out that this pulsar-collapsing dark sector can resolve the core¬ 
cusp problem and will either be excluded or discovered by upcoming direct detection experiments. 
Another implication is a maximum pulsar age curve that increases with distance from the galactic 
center, with a normalization that depends on the couplings and masses of dark sector particles. 
Finally, we use old pulsars outside the galactic center to place bounds on asymmetric Higgs portal 
models. 


I. INTRODUCTION 


Although there is compelling evidence for the exis¬ 
tence of particulate dark matter (DM) from astronom¬ 
ical observations, the precise nature of DM remains elu¬ 
sive. One possibility is that a weakly interacting massive 
particle (WIMP) froze out during the primordial epoch 
with a picobarn annihilation cross-section. However, the 
WIMP paradigm does not explain why the relic abun¬ 
dance of DM is comparable to that of baryonic matter: 

mm- The baryonic relic abundance itself 
was determined by a baryon asymmetry arising at some 
temperature hotter than the baryon-anti-baryon freeze- 
out temperature in the early universe, T > 10 MeV. 
While the exact source of the baryon asymmetry is un¬ 
known, the Sakharov conditions stipulate that particle 
charge (C) and charge-parity (CP) must be violated out 
of thermal equilibrium to generate a particle asymmetry. 
There are numerous models that meet these conditions 
for baryons, including out of equilibrium decays mm and 
the Affleck-Dine mechanism [^. 

Because the source of C and CP violation responsible 
for the baryon asymmetry is unknown, it is plausible, and 
indeed the coincidence Ddm suggests, that the 

dark sector participated in the generation of the baryon 
asymmetry. In the simplest dark asymmetry models (see 
e.g. [IDHE]), DM is charged under a quantum number 
that is composed of or identical to baryon number, and 
the DM abundance is set by the same process that sets 
the baryon abundance. Such models usually imply a DM 
mass range mx = 1 — 15 GeV and collectively fall un¬ 
der the rubrik of asymmetric dark matter (ADM) (see 
Refs. [TMT5] for recent reviews). 

After the dark asymmetry is generated, the final ADM 
relic abundance will depend on having a large annihila¬ 
tion cross-section at freeze-out. ADM freeze-out cross- 
sections must be larger than that of WIMP dark mat¬ 
ter, because if the dark asymmetry provides for most 
of the dark matter relic abundance, ADM annihilation 
cross-sections must exceed the standard picobarn-size 
DM freeze-out cross-section. Otherwise, ADM would 
freeze-out to an overabundance and collapse the uni¬ 
verse. Altogether, these relic abundance requirements 
motivate ADM coupled to light mediators, because a 


picobarn cross-section is difficult to attain using heavy 
(>100 GeV) mediators, while remaining consistent with 
collider bounds muzQ on low mass dark matter (i.e. 
mx = 1 — 15 GeV). However, collider bounds on light 
DM can be evaded if DM annihilates to SM particles 
through a light mediator mixed with a SM boson m- 

In addition to evading collider bounds and matching 
cosmological requirements, light mediator dark matter 
also fits a number of astrophysical anomalies [IH [201126] , 
including the core-cusp problem - some dwarf galaxies 
exhibit unexpectedly cored DM density profiles HZlllH] 
- and the too big to fail problem ~ there is a dearth of 
massive Milky Way satellite galaxies [2^31] . 

More recently, radio surveys of the galactic center have 
discovered a missing pulsar problem [351 [33] . Based on 
the number of progenitor stars which collapse into pul¬ 
sars, around 10 young pulsars were expected to have al¬ 
ready been found within the Milky Way’s central parsec. 
However, so far no young pulsars have been observed. 

In this paper, we demonstrate that the missing pul¬ 
sar problem specifically motivates asymmetric fermionic 
dark matter coupled to the SM through a light Higgs por¬ 
tal mediator, which would collect in and collapse pulsars 
at the galactic center, yet remain consistent with older 
pulsars seen in less DM-dense regions. The structure of 
this paper is as follows: in the remainder of the introduc¬ 
tion, we provide a review of the missing pulsar problem. 
In section |H| we introduce the Higgs portal model, and 
catalogue constraints on it from relic abundance, primor¬ 
dial nucleosynthesis, direct detection, and Higgs invisible 
width measurements. In section IHII we determine what 
Higgs portal parameter space could simultaneously ex¬ 
plain the absence of pulsars at the galactic center, while 
remaining consistent with old pulsars observed near the 
solar position. Section IV gives the expected maximum 
age of pulsars as a function of distance from the galac¬ 
tic center. We present concluding remarks in section [V] 
Appendix|^details estimates of the expected pulsar pop¬ 
ulation at the galactic center. Appendix [B] addresses the 
velocity dependence of the DM-nucleon cross-section. 


^ DM-nucleon couplings for mx = 0.01 — 1 GeV and sub-TeV 
mediators are most constrained by quarkonium—>invisible decays 
at Belle and BES 1181 . 
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A. The missing pulsar problem 

Radio surveys of the galactic center have not found as 
many pulsars as expected [32j . Pulsars are rapidly rotat¬ 
ing, strongly magnetized neutron stars with spin period 
P and approximate spin lifetime Tp ^ P/2P, formed in 
the supernova collapse of (8-20 solar mass) heavy pro¬ 
genitor stars (see [31] for a review). There are about 300 
heavy progenitor stars within a parsec of the galactic cen¬ 
ter [3S], which signify the presence of at least 500 young 
(0.01-100 Myr) pulsars in the same region [3^l33l[36] , 

In addition to these young pulsars, one thousand old 
(^ Gyr) millisecond pulsars (MSPs) are expected to re¬ 
side in the central parsec. Millisecond pulsars form when 
a neutron star in a binary system accretes material and 
angular momentum from its binary companion, thereby 
spinning up to millisecond pulse periods. The dense stel¬ 
lar environment of the galactic center should host many 
binary systems; extrapolating from the population of 
millisecond pulsars observed in star-dense globular clus¬ 
ters (e.g. TerzanS has ^ 100 MSPs [37]), an estimated 
500 — 10"^ millisecond pulsars should abide in the central 
few parsecs [35] . 

In addition to young and millisecond pulsars, there are 
also magnetars, which have a > 100 times stronger mag¬ 
netic field and slower pulses than most pulsars. To date, 
a single magnetar 0.1 parsecs from the galactic center, is 
the only pulsar discovered inside the central 50 parsecs 
[39] [40]. Prior to the measurement of this magnetar’s 
radio pulse dispersion [H], it was reasonable to assume 
that pulsars had not been seen in the galactic center, be¬ 
cause their radio pulses were smeared out by compton 
scattering off a dense intervening column of electrons. 

However, using radio pulses from the magnetar, the 
density of electrons between the GC magnetar and earth 
has been measured. This measurement implies that, with 
the most conservative assumptions, prior radio surveys of 
the galactic center [H] |33] should have found at least 10 
young pulsars and 4 millisecond pulsars [33]. (Note that 
only about a tenth of existing pulsars will have pulses 
that beam towards earth.) A Bayesian analysis of the 
missing young pulsar problem found that the 6.6 GhZ 
galactic center survey [H] puts a 99% upper limit of 200 
young pulsars in the central parsec, comparable to 500 
expected using high mass progenitor data [33]. In Ap¬ 
pendix [^ we provide more detail on methodologies for 
estimating the galactic center pulsar population. There 
are a number of possible explanations for the missing pul¬ 
sar problem - the 300 heavy progenitor stars observed 
inside the central parsec could be a recent anomaly, and 
not indicative of heavy progenitor populations 100 Myr 
ago [32]. Some estimates of the galactic center millisec¬ 
ond pulsar and young pulsar populations assume a typ¬ 
ical stellar initial mass function in the galactic center 
(see Appendix [A|. The actual initial mass function of the 
galactic center may be “top heavy” [45] , resulting in more 
black holes and fewer pulsars produced through core col¬ 
lapse of GG high mass progenitor stars. In this paper. 


we focus on a different, intriguing possibility: asymmetric 
fermionic dark matter coupled to SM particles through a 
Higgs portal may collapse pulsars in the galactic center. 

II. HIGGS PORTAL MEDIATORS FOR PULSAR 
COLLAPSE 

While a number of astrophysical anomalies motivate 
fermionic dark matter coupled to an MeV-scale vector 
or scalar mediator [n mss], in the case of pulsar¬ 
collapsing dark matter (PCDM), scalar mediators and 
asymmetric DM are required. The dark matter must 
be asymmetric, so that it can collect in pulsars with¬ 
out annihilating to SM particles gulls]. If it is coupled 
to a light scalar mediator, it will only have attractive 
self-interaction^ permitting DM self-attractive forces to 
overcome Fermi degeneracy pressure and initiate black 
hole collapse at the center of pulsars. 

Assuming the dark mediator </) is a real scalar, and X 
is a dirac fermion, we add the following terms to the SM 
Lagrangian, 

£ = rsM + iX^X + ^(a</>)2 - acj,{\H\^ - uV2) 

- - uV2) - gD(l)XX - - mxXX 

( 1 ) 

Thus (p and X interact with the visible sector through 
H — (j) mixing [50H52]. Note that we have taken H —>■ 
h + vl\f2 as the Higgs vacuum expectation value (vev) 
after electroweak symmetry breaking and have shifted 
the H — (f> coupling so that p does not get a vev. In the 
low momentum limit, and assuming a,m^ <C v,mh, the 
effective coupling of (p to standard model particles is 

eh = V2av I ml. ( 2 ) 

The coupling of the light mediator to DM, gjy (or 
equivalently ao = g'Dl^'’^), is bounded from below by 
relic abundance and GMB constraints. For thermal 
dark matter, relic abundance requirements put a lower 
bound on the annihilation cross-section of dark matter 
{av) > 10“^®cm^/s, so that the universe does not col¬ 
lapse. If the scalar mediator is lighter than dark matter, 
this implies 



1/2 


\ 10-25cm3/s^ 

1 VGcV/ 

V20/ 


(3) 


^ Note that DM-mediator couplings of the form g£)X'y ^will 
be repulsive for XX XX and XX XX scattering, 
whereas scalar couplings g^XXcj) produce purely attractive self¬ 
interactions. In both these cases the potential is Yukawa \V\ = 
g^exp{—m^r)/Anr. See |49| for a treatment of pseudoscalar and 
axial couplings, which we do not consider here. 
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where mx is the mass of DM, Xf = mx/Ipo, and Tpo 
is the freeze-out temperature of DM [19] . For asymmet¬ 
ric dark matter, requiring that > 0.11 (since the 

relic abundance is assumed to be generated by a parti¬ 
cle asymmetry), the annihilation cross-section should be 
(av) > 4.5 X 10“^®cm^/s for Dirac fermions |53|. 

Because the mediator (f) must decay before big bang 
nucleosynthesis, there is a lower bound on its coupling 
to SM particles. Big bang nucleosynthesis (BBN) takes 
place in the CMB radiation era after about a second. In 
this era, any extra particles that can decay and inject 
energy lead to an enhancement in the nuclear interac¬ 
tion rate. In particular, a light mediator with a mass of 
1 —100 MeV is below the "^He binding energy; therefore it 
would disassociate deuterium and lithium isotopes, while 
leaving "^He relatively undisturbed. Thus there is a con¬ 
straint on any light degree of freedom besides those of 
the SM [54]. To avoid these constraints, we require that 
(p decay before BBN. For < 2 to^, the total width of 
4> is 




epnlm^ 


167rt)2 

2 2?, 



Ai -b 

512ttv^ 



3/2 

0 (to 0 - 2me) 


^ 1/2 + 



(4) 


where the second contribution is for cp —>■ 77 , with Ai = 
— 7 from the W loop, and A^/i = 4/3 is the contribution 
from heavy quark loops [55] . In Figure the dark cyan 
dotted line indicates parameter space where the lifetime 
of p is less than a second. Note that even this constraint 
can be lifted entirely by coupling p to additional light 
dark states (e.g. sterile neutrinos [55]1. 

p — H mixing also contributes to the invisible Higgs 
decay width. The additional contributions to the Higgs 
invisible decay width from Q are 

_ apelmh A ^ / 4m^ 

2 \ ml J S-Kmh y ml 

(5) 



Requiring the Higgs—J-invisible branching ratio be less 
than 40% (from e.g. [^) of the total SM Higgs width 
{Th ^ 4.1 MeV) puts an upper bound on eh- Since 
we only consider m^ m/j, the constraint on the higgs 
invisible decay requires b 0.013. Assuming the contri¬ 
bution from b is negligible, in Figure [^ we indicate the 
constraint on Ch from restricting the Higgs invisible width 
to be less than 40% of the SM expectation [55H5U| . 

DM-nucleon scattering experiments such as LUX m, 
XENONIOO [12], and SuperCDMS [IS] also constrain 
Higgs portal parameter space. A calculation of DM- 
nucleon scattering is given in Appendix [^ in which we 
include the full dependence of the scattering amplitude 
on m^ and the velocity of DM in the nucleon’s rest frame. 
Figure[^shows the current bounds [611163] (solid magenta 
line) and future reach [M] (dashed lines) of DM direct de¬ 
tection experiments. 


III. COLLECTION AND COLLAPSE OF HIGGS 
PORTAL DARK MATTER IN PULSARS 

Another possible method for detecting asymmetric 
Higgs portal dark matter, is to find pulsars imploding 
(or having imploded) in regions where ambient dark mat¬ 
ter is dense and can rapidly collect in pulsars. After 
enough asymmetric dark matter has been captured and 
thermalized in a pulsar, it may form a black hole that 
swallows the pulsar [471 EHl I65U85] . Usually, fermions 
will not form a black hole until N^h ~ {ixipi/mx)^ 
particles have agglomerated, where Nch is the Chan¬ 
drasekhar limit for fermionic matter. The asymmet¬ 
ric Higgs portal model we consider has attractive self¬ 
interactions (from gjjpXX) that counteract Fermi de¬ 
generacy pressure, effectively decreasing Nch- Indeed, 
attractive self-interactions are required for fermionic DM- 
induced pulsar collapse, because the maximum num¬ 
ber of DM particles captured by a neutron star in a 
gigayear, ^ 10 '‘^(GeV/mx), is much smaller than the 
Chandrasekhar number, ^ 10®^(GeV/mx)^. 


A. DM collection 


Dark matter particles are captured in neutron stars at 
a rate given by [HI EH |66l HZ] , 


./6 

^ Px^ 

NbP.vI,c 

V TT ' 

\Vx) 

mx 


I - 


I - exp(-B2) 
52 


ficTux), 

( 6 ) 


where px is the DM density around the neutron star, 
vx ~ 200 km/s is the velocity dispersion of the dark 
matter-neutron star system, Nb = 1.2 x 10 ^^ is the num¬ 
ber of nucleons in a 1-AMq neutron star, Vesc — 0.7 is 
the escape velocity from the surface of a neutron star, 
and ^ = Min[mx/(0.2 GeV),l] is a factor account¬ 

ing for Pauli blocking of DM-nucleon scattering. Be¬ 
cause the dark matter will have a semi-relativistic mo¬ 
mentum when falling through the neutron star’s rest 
frame, the size of its momentum px mx indicates 
that heavy (mx ^ itib) dark matter may not trans¬ 
fer enough momentum to be captured by the neutron 
star. The term in square brackets, with the quantity 
= 6vlgc''^x'mB/vj^imx + mB)^, accounts for this 
diminution in capture of heavier dark matter. Finally, 
the term dependent on the dark matter-nucleon scatter¬ 
ing cross-section, a^x, is given by f{anx) = crsat(l - 
exp(—cTiix/crsat))- When the dark matter-nucleon scat¬ 
tering cross-section is small, anx 10 “^® cm^, this 
function returns anX, but as (t„x gets larger, the geo¬ 
metric cross-section per nucleon in the neutron star sat¬ 
urates. Note that for lighter dark matter, this satura¬ 
tion cross-section will depend on Pauli blocking [48l HZ] j 
<Xsat = R^s/0.45Vb^, where Rxs ~ 10 km is the neu¬ 
tron star radius. We detail our calculation of the dark 
matter-nucleon scattering cross-section cr„x for a given 
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px = 80 or 0j2 GeV/cm^, tns = 0.01 or_l. Gyr, v = 150 or 200 km/s, 
Tns= 10^ K, Mns= 1.4 M-, Rns= 10 km 




Figure 1. Direct detection, cosmological, and collider bounds 
on Higgs portal dark matter are shown as a function of th = 
y/2avlm\, along with viable galactic center pulsar-collapsing 
parameter space, shaded grey. Space above the solid black, 
blue and pink lines are excluded by Gyr old pulsars outside 
the GC, constraints on the Higgs invisible decay width, and 
direct detection, respectively. Constraints from BBN exclude 
space below the dotted cyan line, but these can be removed 
by adding decay modes for cj), Ref. [Sn|. The DM density 
(px), pulsar lifetime (tns), velocity dispersion (u), and pulsar 
temperature {Tns), mass {Mns), and radius (Rns) used to 
plot the solid and dotted-dashed black lines are given at the 
top of the plot. Note that parameter space along the dotted- 
dashed (solid) black line would collapse pulsars older than ~ 
10^ (~ 10®) years inside the Milky Way’s central 500 parsecs, 
while permitting up to ~ 10^® (~ 10®) year old pulsars outside 
of the galactic center. 


set of Higgs portal parameters {eh,rnx,'m^,aD) in Ap¬ 
pendix 


B. DM state at collapse 


To determine the critical number of dark fermions that 
will initiate black hole collapse in the neutron star inte¬ 
rior, we use the Virial equation for a test particle at the 
edge of a sphere of fermions thermalized at the center of 
the neutron star. 


—2Ek + 


7 ^ 


+ Vvuk = 0 , 


(7) 


where the first term is the virialized kinetic energy, the 
second term is the virialized gravitational potential from 
baryons in the neutron star, and the third term is the 
virialized form of the Yukawa potential. It is reason¬ 
able to assume all DM particles are thermalized, be¬ 
cause results from [Eli US] indicate that for the param¬ 
eter space under consideration, dark matter thermalizes 
inside the neutron star on timescales < 10® yrs. Nx 
is the number of dark matter particles collected, and 
Vyuk = Y.rj -b m^), with rj being the 

inter-particle distance. Once Vyuk ^ ‘^Ek, the dark mat¬ 
ter will collapse into a black hole. In these calculations, 
we can neglect the dark matter self-gravity as being small 
(in what follows, see [SI] for more details). In Eq. ([^, 
we have defined y = which is the expo¬ 

nent of the Yukawa potential if only the nearest fermions 
to the test particle contribute to its potential; note that 
X = 1.6r/iV^^^ is the nearest-neighbor inter-fermion dis¬ 
tance for Nx fermions evenly distributed in a sphere of 
radius r. 

Strongly-screened. If y 3> 2, only the nearest-neighbor 
particles will matter in determining the Yukawa poten¬ 
tial. We will call this the “strongly-screened” limit, the 
limit in which all but the nearest-neighbor fermions can 
be neglected. The strongly-screened virialized Yukawa 
potential is 

{mtj,e~y/y + , (8) 

where we assume eight nearest neighbor particles. 

Coulombic. If y <C 2, on the other hand, the expo¬ 
nential piece of the Yukawa potential approaches unity, 
and the Yukawa potential becomes Coulombic. In other 
words, the Yukawa potential (oc will not be 

suppressed inside radius l/m^. The number density of 
dark matter fermions at the center of the star is given 
by 1/x®. Hence, the number of fermions contribut¬ 
ing to the Coulomb-like potential inside radius l/m^, is 
Nco = 47r/3m^a:®, which gives a virialized potential term, 

^Yuk’' = iaoNcom^, = ^TTaom^/y^. (9) 

Degenerate. Before the onset of collapse, the dark mat¬ 
ter fermions collected inside radius r will be degenerate 
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if more than 

Ndeg = 5 X 1027(r/cm)3(mx/GeV)3/2 (10) 

have collected, assuming an ambient temperature of T = 
10"^ K. In this case the kinetic energy is given by 

» = (9^iVx/4)"/V2mxr" = {ZT^y^m\l2mxy^ 

( 11 ) 

and the pre-collapse radius of the dark matter, deter¬ 
mined by solving Eq. Q with the last term omitted is 

rthMg = 2.4 X 10"^7Vj/®(GeV/mx)^/^ cm. (12) 

Non-degenerate. If less than Ndeg particles have 
collected before collapse, the kinetic energy is simply 
j^non-deg ^ and the pre-collapse thermal radius 

is 

^th,nondeg — 80(GeV/mx)^^^ cm, (13) 

again assuming a neutron star temperature of 10^ K. 

C. Critical DM number for collapse 

There are four states from which the collapse of DM 
can begin: the dark matter may be either degenerate or 
non-degenerate and either strongly-screened or Coulom- 
bic. We can determine which state precedes collapse 
as a function of model parameters (aDjirix^rn^) using 
Eq. 0. 

Degenerate. If more than N^eg particles have collected, 
the dark matter will be degenerate. Substituting 
into Eq. and dropping the baryonic gravity term 
(which should be negligible after collapse begins), we can 
determine what parameter space initiates Goulombic col¬ 
lapse by requiring y <2 in 

(37r2)2/3^y2mxj/^ = AiraDm^/y^, (14) 

from which we find that collapse will begin from a 
Goulombic (strongly-screened) state when a > m^/mx 
{a < m^/mx). If collapse begins from a Goulombic, de¬ 
generate state, the number of particles required for col¬ 
lapse can be found by solving Eq. HI for Ncou, 

> 10^^a-®(m,^/MeV)^2(GeV/mx)®. (15) 

If on the other hand, collapse begins from a strongly- 
screened degenerate state, Eq. 0 must be solved nu¬ 
merically with ■ Specifically, we solve Eq. 0 for 

y, using successively larger values of Ncoii, until no so¬ 
lution can be obtained where y > 2, indicating collapse 
has begun. More details of this numerical solution can 
be found in [82] . 

Non-degenerate. If less than Ndeg particles have col¬ 
lected, dark matter will not be degenerate, and the ki¬ 
netic energy will be 3kBT/2. For all parameter space 



mx (GeV) 




Figure 2. Constraints from SuperCDMS and LUX (pink) 
along with bounds from old pulsars (black) on asymmet¬ 
ric Higgs portal dark matter are shown as a function of 
= y/^avfm^, for a number of values of an and m^. The 
DM density {px), pulsar lifetime (tjvs), velocity dispersion 
(v), and pulsar temperature (Tjvs), mass (Mjvs), and radius 
(Rns) used to plot the black lines are given at the top of the 
plot. Note that parameter space along each black line would 
collapse pulsars older than ~ 10® years inside the Milky Way’s 
central 500 parsecs, while permitting ~ 10® year old pulsars 
outside of the galactic center. 


of interest (mx ^0.1 GeV), the number of particles re- 
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quired to begin collapse in a Coulombic, non-degenerate 
state exceeds Ndeg [22] ■ Hence, we only need to consider 
collapse from a strongly-screened non-degenerate state. 
To determine the critical number of particles required, we 
insert Q and solve numerically for Ncoii- 

As collapse proceeds, and y (along with r) shrinks, the 
fermions will eventually become degenerate. The value of 
y at which the fermions become degenerate can be found 
by substituting rth,deg into y, 


ydeg = 2 O(m 0 /MeV)(GeV/mx)^/^. (16) 


If ydeg < 2, we can substitute Eq. (16 1 into Eq. (14), 
and find that collapse continues through a Coulombic, 
degenerate state so long as 


a > 0.02(TO,^/MeV)2(GeV/mx)^/^. (17) 

If Vdeg > 2, then substituting Eq. into Eq. (drop¬ 
ping the baryonic term), we find that collapse proceeds 
so long as 


a > 10 ®e*'‘^'s(MeV/m0)(l-f 1/ydeg) ^ (18) 


After dark matter collapses into a black hole, in order 
to destroy the neutron star, it must accrete surrounding 
baryons faster than it radiates via the Hawking mecha¬ 
nism [87]. This condition is fulfilled if, 

(19) 


where the first term is the Hawking radiation rate, the 
second term is the Bondi accretion rate, and Cg 0.3 
is the sound speed of baryons in the neutron star. If 
the black hole grows, it will destroy the neutron star in 
tdes ~ 2000(GeV/mx)(1047Ncou) yrs. 

In Figure [2 and we display parameter space which 
is excluded by the existence of 10 ® year old pulsars 
outside the galactic center (where the DM density is 
px 0.2 GeV/cm®). Note that the left side of these 
bounds terminate where Ncou. > Cx^nSi while the right 
side terminates when the condition Eq. (19) fails, mean¬ 
ing the black holes formed by collapsing dark matter are 
too small too grow and instead evaporate. In Figure [l] 
we show model space which would collapse 10 ® — 10 ^ year 
old pulsars in the galactic center, while allowing for old 
pulsars outside the galactic center. Note that we assume 
a galactic center DM density {px 80 GeV/cm®) deter¬ 
mined by DM self-interactions, as we discuss in the next 
section. 


IV. DARK MATTER SELF-INTERACTIONS 
AND PULSAR AGE CURVES 

To determine which Higgs portal PCDM models pre¬ 
cipitate pulsar destruction in the galactic center, we need 
to know the mass distribution of DM in the Milky Way. 


It is commonly assumed that DM is cold and collisionless 
(GGDM), and simulations of CCDM structure formation 
yield halo densities that follow the Navarro-Frenk-White 
(NFW) profile: p{r) = pQ/{r/Rs) (1-\-r/Rs)~'^, where 
Po and Rs are characteristic density and scale radius re¬ 
spectively, and r is the distance from the galactic cen¬ 
ter. However, the CCDM paradigm has been challenged 
recently by some observations of small scale structure. 
These include dwarf halo prohles which appear under- 
dense at their centers |27l [28], too few small satellite 
galaxies near the Milky Way, and no large satellite galax¬ 
ies p5H5T] . 

These small scale structure anomalies can be addressed 
with the addition of sizable DM self-interactions asi 
10“®"^ cm®/GeV. However, this seems disfavored by self¬ 
interaction bounds from “bullet” collisions in galactic 
clusters jSS] [89| and the observed ellipticity of galac¬ 
tic cluster cores [90]. Typically these require DM self¬ 
interactions smaller than asij'mx 10“®"^ cm®/GeV. 
If, on the other hand, DM self-interactions are velocity 
dependent through coupling to a light mediator (m^ = 
1 — 100 MeV), bounds from large scale observations can 
be evaded mUDHlSI. This is because the cross-section 
can be large in, for example dwarf galaxies, where the 
characteristic velocity is 50 km/s, while the mediator is 
more off-shell and interactions are weaker for DM veloc¬ 
ities ^ 1000 km/s, in galactic clusters. Much of the pa¬ 
rameter space shaded grey in Figure will have velocity- 
dependent self-interactions suitable for resolving small 
scale anomalies, while remaining consistent with large 
scale bounds [ 2 ^ |35] . 

Some simulations indicate that DM with signihcant 
self-interactions a si ^ 10“®"^ cm®/GeV, will cause the 
Milky Way’s DM density profile to be flat (~ GeV/cm®) 
within the central 2 kiloparsecs pTI193j . However, a re¬ 
cent study which included the gravitational potential of 
baryons showed that SIDM in Milky Way size halos will 
create a 500 parsec DM core, with a larger constant cen¬ 
tral density 7 80 GeV/cm®) (Mj. For SIDM, the den¬ 
sity of the central core will depend mostly on central the 
baryonic density, p(r) = po 6 xp[/i(r)], where h{r) satisfies 
Jean’s equation and is given by 


h{r) 


-27rGpo7 


y 



( 20 ) 


In the Milky Way, po = 80 GeV/cm® is the central DM 
density, y = r/(r + tq) parameterizes distance from the 
galactic center, tq = 2.7 kpc, ctq = 160 km/s is a constant 
velocity dispersion. We show this DM density profile 
in the upper left inset of Figure along with the DM 
velocity dispersion we assumed, which was derived by 
fitting observed star velocities in the Milky Way |84l 
In Figure we use the DM density profile in Eq. (20) 
to predict the maximum pulsar age as a function of dis¬ 
tance from the galactic center. The solid (dotted-dashed) 
line shows the maximum age curve for asymmetric Higgs 
portal models that predict > 10 ® (> 10 ^) year old pulsars 
will collapse in the galactic center. The radial distance 
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Radius from galactic center (kpc) 


Figure 3. The maximum age of an Rns = 10 km, Mns = 
1.4M0, Tns = 10'* K pulsar is plotted as a function of dis¬ 
tance from the galactic center. The solid and dotted-dashed 
maximum age curves result from parameter space shown in 
Figure also indicated in Figure plots with solid and 
dotted-dashed curves. For example, the solid curve is consis¬ 
tent with light dark matter {mx = 6 GeV) strongly coupled 
(od = 0.1) to a light scalar {m^ = 15 MeV) with a small 
mixing {eh — 10“®) with the Higgs boson; this parameter 
space lies along the solid black curve in the top panel of Fig¬ 
ure Pulsar data points are overlaid; “characteristic” ages 
(Tp < P/2P) and radial distances were taken from the ATNF 
pulsar catalog [96]. Note that characteristic ages should be 
treated as upper bounds on pulsar lifetimes - actual pulsar 
ages can be significantly shorter, especially if Pq ~ P where Po 
is the initial pulse period. (Indeed, many pulsar characteristic 
ages above the dotted-dashed curve exceed the age of the uni¬ 
verse.) Pulsar data points shown in dark blue were required 
to have longer pulse periods P > 0.5 s and to be solitary (not 
in a binary); these conditions may indicate better alignment 
between characteristic and actual age [34]. Inset are the dark 
matter density (top) and velocity dispersion (bottom) as a 
function of distance from the galactic center. 


from the galactic center and the characteristic radii of 
pulsars in the Milky Way, shown with blue and pink data 
points, are taken from the ATNF pulsar catalogue [96]. 
It should be pointed out that the characteristic age of 
pulsars Tp < Pl2P, often vastly overestimates the age of 
the pulsar. Indeed, it has become clear that this over¬ 
estimate may be systematically large for all millisecond 
pulsars with Tp > 10® years [331117]. This effect results 
from an initial pulse period, Pq, that is not much smaller 
than the current pulse period P. In this case the full 
expression for the pulsar age, ~ P/2P — Po/2P, is 
required, but Pq is unknown. Note in Figure[^that many 
of the pulsar “characteristic” ages exceed the age of the 
universe. 

However, millisecond pulsars that formed in a binary 
system with a white dwarf (WD) provide an alternative 
way to check the age of the pulsar. The binary white 
dwarf’s temperature and mass can be fit to WD cooling 
curves to independently determine the age of the binary 


system. One example is pulsar J1738-I-0333, which has a 
characteristic age of 4 Gyr, and whose white dwarf has an 
apparent age of 0.5 — 5 Gyr [15]. Therefore, for the pur¬ 
pose of bounding DM parameters with pulsar ages, it is 
most conservative to assume pulsars outside the galactic 
center have reached ages of about a gigayear. 


V. CONCLUSIONS 

Higgs portal DM-induced pulsar collapse could ex¬ 
plain the apparent paucity of galactic center pulsars. 
Fermionic asymmetric dark matter coupled to the SM 
through a Higgs portal can collapse pulsars older than 
~ 10® years inside the galactic center, while not collaps¬ 
ing 10® year old pulsars near the solar position. This is 
because the expected DM density in the galactic center 
{px = 80 GeV/cm®) is larger than in the rest of the halo 
{px = 0.2 GeV/cm®), leading to increased DM capture 
in GC pulsars. 

Based on the number of GG pulsar progenitor stars, 
GG radio surveys should have already found 0(10) pul¬ 
sars in the central parsec. However, none have been ob¬ 
served, and it is unlikely that this is solely the result 
of measurement limitations |32j . although it could result 
from an overestimated GC pulsar population. We have 
showed that for DM masses mx = 0.1 — lOOGeV, with 
Higgs portal mediator masses 5 — 20MeV, DM-mediator 
couplings an = 0.001 — 0.1, and mediator-Higgs mixings 
eh = 10“® — 10“^ asymmetric Higgs portal dark matter 
provides an explanation for the absence of GC pulsars: 
older GC pulsars collapse into black holes after enough 
DM is captured to form a black hole that will grow in 
pulsar interiors. As pulsars are discovered in the galactic 
center by future radio surveys and radio telescopes, (e.g. 
FAST [99] and the Square Kilometer Array [100] ) PCDM 
would manifest as a maximum pulsar age that increases 
with distance from the GC. 

Applied to old pulsars seen outside the GC, our results 
set bounds on fermionic, asymmetric Higgs portal dark 
matter that are often more stringent than those set by di¬ 
rect detection experiments. However, while much of the 
Higgs portal PCDM parameter space is inaccessible to 
present terrestrial direct detection, we have shown that 
next-generation direct detection experiments will either 
find or exclude it. In addition, through its coupling to 
a light mediator, Higgs portal PCDM naturally fits one 
explanation for the core-cusp phenomenon: light media¬ 
tor DM self-scattering is resonantly enhanced at smaller 
momenta in dwarf galaxies, and is diminished at larger 
momenta in spiral galaxies and galactic clusters. 

Finally, we note that a recent proposal of Fuller and 
Ott [101] (advanced with all reserve) has linked PCDM 
to fast radio bursts [102H106] . We leave this and other 
complementary probes of PCDM to future work. 
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Appendix A: Notes on the galactic center pulsar 
population 


This appendix details methods for estimating the num¬ 
ber of pulsars, both young and millisecond, at the galactic 
center. 

Young pulsars in the central parsec. A simple way to 
estimate the expected number of young pulsars in a re¬ 
gion is to count the number of 8-20 solar mass, high mass 
stars, Nhrns- It is expected that the majority of these 
will expire in core collapse supernovae and form pulsars 
[82l [36] . More massive stars have larger hydrogen cores 
which burn faster, and so as a general rule the lumi¬ 
nosity of stars scales as L ~ Normalizing to a 

sol-type star, the lifetime of a high mass star (t ~ M/L) 
is given by thms = yrs. Assuming that 

the abundance of each star type stays constant in the re¬ 
gion of interest, this implies there are ~ hNhms/'i young 
pulsars, because the lifetime of the high mass stars is 6 
Myr, and a typical lifetime for a young pulsar is 10 Myr. 

It follows that the inner parsec of the galactic center 
which hosts ~300 high mass stars should also host ^500 
young pulsars, ^50 of which would beam towards earth, 
^25 of which should have already been detected given the 
recently measured radio scattering dispersion (courtesy 
of a bright magnetar) [52] ■ However, this calculation is 
only valid if high mass stars in the galactic center were 
as abundant historically as they are now; high mass stars 
may have been less abundant 20 million years ago [35]. 
The next generation of surveys should permit imaging 
of 1-2 solar mass stars m. which as the byproducts of 
high mass stars, will provide a better indication of the 
historical abundance of high mass progenitor stars. 

Older, millisecond pulsars. Millisecond or recycled 
pulsars form when a neutron star spins up in a binary 
system by accreting gas from its companion [33]. Most 
millisecond pulsars have been found inside globular clus¬ 
ters, which have stellar densities as large as 1 — 10^ solar 
masses per cubic parsec, comparable to a stellar density 
of lO^Mo/pc^ and lO^Mo/pc^ in the central 100 par¬ 
secs and central parsec of the milky way, respectively. 
One way to estimate the central parsec millisecond pul¬ 
sar population is to scale up the millisecond population 
in globular clusters (about 10-50 millisecond pulsars per 


cluster) [108] . Millisecond pulsar formation is expected 
to be greater in the central parsec, because of its higher 
stellar density and higher escape velocity (400 vs. 50 
km/s), allowing it to retain a larger fraction of young neu¬ 
tron stars (with an average birth velocity of 400 km/s). 
Given these factors, ~1000 milliscond pulsars are ex¬ 
pected in the central parsec [551 110511109] . 

We can also consider a simpler, more conservative 
method for estimating the millisecond pulsar population 
of the central parsec. Before a neutron star has finished 
accreting gas from its companion, it will necessarily be 
in a low mass x-ray binary. As gas falls into the compact 
neutron star, it often emits x-rays. Hence, we should 
expect the observation of low mass x-ray binaries to cor¬ 
relate with the number of millisecond pulsars. Indeed, 
studies have shown that both of these correlate with the 
stellar encounter rate in globular clusters, Tc pcr^v~^, 
where these variables are respectively the density, ra¬ 
dius, and linear velocity dispersion of a globular clus¬ 
ter [nniiiii]. Roughly 10 — 20 times more millisecond 
pulsars than low mass x-ray binaries have been found in 
globular clustars m- Extrapolating from the four x-ray 
binaries found in the central parsec mm, one would 
expect ~40-80 detectable millisecond pulsars in the cen¬ 
tral parsec, of which (accounting for increased radio pulse 
dispersion at the galactic center) ~5-10 would have been 
detected by the 14 Ghz survey [321 112 ■ 


Appendix B: Onx velocity dependence for Higgs 
portal DM 


The matrix element for dark matter t-channel scatter¬ 
ing off nuclei is 


iM = gaf-N 


[u{pa)u{pi)][u{p3)u{p2)\ 
t — -I- ie 


(Bl) 


where and p 2,4 are the dark matter and nucleon ini¬ 
tial and final 4-momenta, respectively. The prefactor 
Eat — 3 X 10“^ X Ch parameterizes the coupling of </> to 
nuclei, largely determined by the gluon hadronic matrix 

element m- 

In the center momentum frame of the interaction, the 
four momenta are given by 

Pi — (Ei,0,0,Pcm) 

P2 = (^ 2 , 0 , 0 , -Pcm) 

P3 = (Es, -PcmSin 0, 0, -pcmCOS 0) 

P4 = (E4,PcmSin 0,O,PcmCOS 0). (B2) 


To approximate the cross-section at a direct detection 
experiment or in a neutron star, we integrate the square 
of the amplitude, 

l\Mf = 

2 2 4(2m^ -bpg^(l - cos 0)){2ml + pl^jl - cos 0)) 

(-2pL( 1 - cos 0) - mlY 

(B3) 
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over a range of nuclear recoil energies i?recoii, specifi¬ 
cally ranging over incoming DM momenta O.OOlyumjf to 
'yvmx, where 7 = 1/Vl — and u is a typical DM ve¬ 
locity in the nucleon’s rest frame {v ~ 10“^ c and 0.7 c for 
DD experiments and pulsars, respectively). Here we de¬ 
fine the momentum transfer Q = ~ which 

can be related to the recoil energy with Pcm(l ~ cos 6) ^ 


Some recent work has shown that Dirac dark matter 
with a strong (a_D > 0.1) coupling to a light scalar me¬ 


diator, may form bound states at freeze-out [53 I115j . 
This bound state dark matter can have a cross-section 
enhanced by the square of the number of bound 
fermions, as long as the momentum transfer of the inter¬ 
action is less than the binding energy Q ^ BEq. For 
the purposes of this study, we assume no bound state en¬ 
hancement to the cross-section - this assumption is cer¬ 
tainly valid for the capture of DM in neutron stars where 
Q ~ mx BEq. However, there may be some enhance¬ 
ment to direct detection capture when ^ Am^/mx, 
which would shift direct detection bounds on in Fig¬ 
ures |T] and 121 
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